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THE YOUNG PLANET-MASS OBJECT 2M1207b: 
A COOL, CLOUDY, AND METHANE-POOR ATMOSPHERE 
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ABSTRACT 

The properties of 2M1207b, a young (~ 8 Myr) planet-mass companion, have lacked a satisfactory 
explanation for some time. The combination of low luminosity, red near-IR colors, and L-type near- 
IR spectrum (previously consistent with T e g ~ 1600 K) implies an abnormally small radius. Early 
explanations for the apparent underluminosity of 2M1207b invoked an edge-on disk or the remnant 
of a recent protoplanetary collision. The discovery of a second planet-mass object (HR8799b) with 
similar luminosity and colors as 2M1207b indicates that a third explanation, one of a purely atmo- 
spheric nature, is more likely. By including clouds, non-equilibrium chemistry, and low-gravity, an 
atmosphere with effective temperature consistent with evolution cooling-track predictions is revealed. 
Consequently 2M1207b, and others like it, requires no new physics to explain nor do they belong to 
a new class of objects. Instead they most likely represent the natural extension of cloudy substellar 
atmospheres down to low T e ff and log(g). If this atmosphere only explanation for 2M1207b is correct, 
then very young planet-mass objects with near-IR spectra similar to field T dwarfs may be rare. 
Subject headings: planetary systems — stars: atmospheres — stars: low-mass, brown dwarfs, stars: 
individual (2MASSW J1207334-393254) 



1. INTRODUCTION 

2M1207b is a planet-mass companion orbiting the 
young brown dwarf 2M1207A at a projec ted separation 
of - 46 AU dChauvin et alJ l200l I2005L As a mem- 
berof the ~ 8 Myr old TW Hydra association (|Gizis I 
2002), 2M1207b is one of the youngest companions of its 
mass (~ 2 - 5 Mj up ) ever imaged. Regardless of how 
the system formed, 2M1207b provides a unique look into 
the atmospheric properties of giant planets and brown 
dwarfs in the very early stages of their evolution. 

Soon after 2M1207b was discovered, it was recog- 
nized that its low luminosity and red near-IR colors 
were seemingly incon sistent. With m^ = 16.93 ±0.11 
(|Chauvin et al.l l2004f). a If- band bolometric correction 
of 3.25 ± 0.14 dMamaie k 2005), and a distance of 52.4 
±1.1 pc (|Ducourant et al.1 12008). 2M1207b has a lumi- 
nosity of logL/L = -4.7 3 ± 0.12. With this l uminos- 
ity and an age of 8tj Myr (jChauvin et al.ll2004D . brown 
dwarf cooling tracks (iBaraffe et al.1 2003) predict 4 T c g = 
936 - 1090 K, log(g) = 3.5 - 3.8 (cgs units), radius = 
1.34 - 1.43 l?j U p, and a mass of 2.3 - 4.8 M Jup . The red 
near-IR colors (e.g., H — K = 1.16 ± 0.24; Chauvin et 
al. 2004) and near-IR spectra of 2M1207b, on the other 
hand, are indicative of a mid to lat e L spectral type and 
effective temperatur e of - 1600K (jMohantv et all 120071; 
iPatience et al.l l2010t ). Such a combination of high ef- 
fective temperature and low luminosity would require a 
radius of ~ 0.7-Rj up , about a factor of two below evo- 
lution model predictions. Under the assumption that 

1 Lowell Observatory, 1400 W. Mars Hill Rd., Flagstaff, AZ 
86001 Email: barman01owell.edu 

2 Lawrence Livermore National Laboratory, 7000 East Avenue, 
Livermore, CA 94550, USA 

3 National Research Council Canada, Herzberg Institute of As- 
trophysics, 5071 West Saanich Road, Victoria, BC V9E 2E7, 
Canada 

4 The ranges in predicted values come from the luminosity and 
age uncertainties. 



the higher T c g is correct, two possible explanations for 
the apparent underluminosity of 2M1207b have been put 
forward. An edge-on disk, albeit with unusual charac- 
teristics, could provide the necessary 2.5 mags of gray 
extincti on to accommodate a larger and more reasonable 
radius (|Mohantv et al.ll2~007j ). A more provocative sug- 
gestion requires 2M1207b to be the hot afterglow of a 
recent protoplanetary c ollision allowing it to be simulta- 
neously hot and small (Mama iek fc Mevedl2007l) . 

Both of these early explanations for the observed 
properties of 2M1207b hinge upon the effective tem- 
perature being as high as 1600 K. Earlier estimates of 
T e e were likely lead astray by the lack of color-T c ff re- 
lationships properly calibrated for young, planet-mass 
objects and the lack of model atmospheres spanning 
a broad enough range of cloud and chemical parame- 
ters to encompas s objects like 2M1 207b. The discov- 
ery of HR8799b (jMarois et al.l I2Q08T) . a second planet- 
mass object with similar near-IR colors and luminosity 
as 2M1207b, motivates this new model atmosphere study 
of 2M1207b, as it is highly unlikely that either the edge- 
on disk or recent-collision explanation applies to both 
objects. 

In this Letter, an atmosphere-only explanation for the 
observed properties of 2M1207b is presented. A combi- 
nation of clouds of modest thickness and non-equilibrium 
CO/CH4 ratio is shown to simultaneously reproduce 
both the observed photometric and spectroscopic prop- 
erties of 2M1207b, with bulk properties consistent with 
evolution model predictions. Such a n explanation was 
touched upon in several recent papers jCurrie et alJ201ll : 
ISkemer et al1l2011l:lBa"rman et alj|2011l) . but here the at- 
mosphere of 2M1207b is explored in more detail. 

2. CLOUDS 

The properties of brown dwarfs and giant planets are 
known to be influenced by atmospheric cloud opacity and 
it is well established that clouds play a central role in 
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the t ransition from spectral types L to T (|Allard et al.l 
2001). Early work on brown dwarf atmospheres ap- 
proached cloud modeling phenomenologically, parame- 
terizing the problem with a n emphasis on vertical mixing 
(jAckerman fe Marlevll200TI ). Cloudy and cloud- free lim- 
its provide useful insight into the expected spectroscopic 
and photometric trends but often fail, unsurprisingly, to 
match individual brown dwarfs, especi ally in the L-to-T 
transition region (|Burrows 

Figure [1] compares field brown dwarfs to 2M1207b 
and the HR8799 planets in a color-magnitude diagram 
(CMD). 2M1207b is located between the cloudy and 
cloud-free limits and, consequently, one should not ex- 
pect either of these limiting cases to be appropriate when 
modeling its photometric and spectroscopic properties. 
Figure Q] shows the path brown dwarfs or giant planets 
of various effective temperatures would follow in a near- 
IR CMD if the vertical thickness of clouds is allowed 
to continuously increase from zero (cloud-free) to well 
above the photosphere (pure equilibrium clouds). The 
radius for these tracks, 1.4 i?j up , was specifically chosen 
to match the radius predicted for 2M1207b as discussed 
above; however, the paths traced by these tracks are oth- 
erwise independent of evolutionary models. 2M1207b is 
intersected by the T e s= 1000 K cloud-track, which is 
the expected value from evolution models, demonstrating 
that low-temperature cloudy atmospheres can achieve 
very red near-IR colors, even with clouds significantly 
thinner than the extreme limits. 

3. NON-LOCAL EQUILIBRIUM CHEMISTRY 

If the methane-rich atmospheres of mid to late T 
dwarfs were in a pure chemical equilibrium state, CO 
mole fractions would be too small to have a major im- 
pact on their spectra. However, CO has been detected 
in many T dwarfs, suggesting that their atmospheres are 
out of equilibrium poll et al.l ll997t ISaumon et aIll2000L 
120061 iGeballe et aT1l2009f ). The most likely mechanism 
for CO enhancement is vertical mixing from deep layers, 
where the temperatures and pressures are higher and CO 
is in ready supply. At photospheric depths and above, 
the chemical timescale (r c hem) to reestablish an equilib- 
rium CO/CH4 ratio becomes far greater than the mixing 
timescale (r m i x ), thereby allowing larger CO mole frac- 
tions to exist at otherwise methane-dominated pressures. 
Through the same mixing process, N2/NH3 can also de- 
part from local chemical equilibrium (LCE; Saumon et al. 
2006; Hubeny & Burrows 2007) . The standard non-LCE 
model quenches the CO and CH4 mole fractions at the 
atmospheric pressur e (P a ) where T m - lx = r c h m, with T m i x 
computed following ISmithl ()1998l ). Below the quenching 
depth (P > P q ) the mole fractions for CO and CH4 are 
in chemical equilibrium while above {P < P q ) they are 
set to the values at P q . 

The mole fractions for CO and CH 4 , N(CO) and 
N(CH4), at and above P q are very sensitive to the un- 
derlying temperature-pressure (T-P) profile and, thus, 
are se nsitive to gravity, cloud opacity, and metal- 
licitv dHubeny fe Burrowi l2007t iFortnev et all 120081 : 
iBarman et al.ll201lD . Certain combinations of low grav- 
ity and clouds can result in P q below the CO/CH4 equi- 
librium chemistry crossing point (P eq )- When P q is suf- 
ficiently deep and P q > P eqi N(CO) can be quenched 
at its maximum value while N(CH 4 ) is quenched near 
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Fig. 1. — Absolute _ff-band magnitude vs. H-K near-IR color- 
magnitude d iagram for field brown dwarfs ( Lcggctt et al] 120021 ; 
IKnapp et aT1l200ll) . The planets 2M 1207b and HR8799bcd are 
indicated by symbols with 1-cr error bars (see legend). Dotted 
lines show color-magnitude tracks for chemical equilibrium mod- 
els with radius = 1.4 i?j up , log(g) = 3.5, mean particle size equal 
to 5 fj,m and T e fj equal to 700 - 1200 K in steps of 100 K, from 
bottom to top. Cloud thickness increases from left to right. The 
arrows indicate the approximate locations of cloud free (left) and 
extremely cloudy (right) models, with arrow-direction pointing to- 
ward decreasing T c g . The "transition" region, between cloudy and 
cloud-free atmospheres covers a broad range in the color-magnitude 
diagram, beyond what is currently occupied by field brown dwarfs. 



its minimum. When this situation occurs, the non- 
LCE CO/CH4 ratio becomes fairly insensitive to the 
mixing timescale in the radiative zone (determined by 
the adopted coefficient of eddy diffusion, K zz ). This 
situation is similar to the N2/NH3 chemistry where 
the NH3 mole fractions can also be nearly indepen- 
dent of K zz , even in high-gravity cloud-free atm ospheres 
(jSaumon et al.ll2006t iHubenv fe Burrowsll2007ft . 

This atmospheric behavior is highly significant to 
2M1207b as it allows for a photosphere with much higher 
N(CO) and much lower N(CH 4 ) mole fractions at the low 
T e g predicted by evolution models. Previous studies, fo- 
cused primarily on field brown dwarfs, present far less 
severe situations with non-LCE only altering spectra at 
A > 4 /^m wher e the strongest absorptio n bands of CO 
and NH 3 occur (|Hubenv fc Burro ws 2007). At lower sur- 
face gravities, however, strong non-LCE eff ects have been 
shown to extend well in to the near-IR (jFortnev et al.l 
l2lMlBarman et al.ll20ll . 

4. MODEL COMPARISONS 

2M1207b has been observed extensively from the 
ground and space, w ith photometric covera ge between 
0.9 a nd ~ 9 urn (Chauvin et al.l 120041: ISong et al.l 
l200?l iMohantv et all 120071 : iSkemer et al.l 1201 lh . High 
signal-to-noise J, H , and K spectra are also available 
(|Patience et al.ll2010[) . With the goal of finding a model 
atmosphere that agrees with these observations and that 
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has T e ff and log(g) in the range predicted by evolution 
models, a sequence of atmosphere models was computed 
covering T cS = 900 - 1200 K and log(g) from 3.0 to 4.5 
(cgs units). The same inte rmediate cloud (ICM ) and 
non-LCE prescriptions from I Barman et al.1 (|2011l ) were 
used. Synthetic photometry was generated by convolving 
synthetic spectra with filter response curves. J, H, and 
K spectra were produced by convolving the model spec- 
tra with a Gaussian filter matching the observed spectral 
resolution, then interpolated onto the same wavelength 
points. The best-fit was determined by least-squares 
minimization. 

Figure [5] illustrates the basic structure of the model 
that best fits the data (T eff = 1000 K and log(g)= 4.0). 
The atmospheric cloud, composed mostly of Fe and Mg- 
Si grains, has a base at ~ 3 bar and extends upward 
before dropping off in number-density at ~ 1 bar. De- 
spite the rapid drop in number-density, the cloud extends 
across the photospheric depths. Also, P q (~ 3 bar) is well 
below the N(CO) e , = N(CH 4 ) eg point (P - 0.3 bar), 
with the CO mole fractions set to the maximum value 
and CH4 is close to its minimum value. 

While the model cloud and non-LCE properties are de- 
termined by free-parameters, they are likely supported 
by low gravity and efficient vertical mixing. In the con- 
vection zone T m i x cx Hp/V c , where V c is the convective ve- 
locity and Hp is the local pressure scale-height. With K zz 
oc i?p/r m i x oc V c Hp, K zz increases with decreasing grav- 
ity in the convection zone as velocity and scale-height 
increase. The radiative-convective boundary also shifts 
toward the photosphere as surface gravity decreases, fur- 
ther suggesting that vertical mixing in the radiative zone 
near this boundary will also be enhanced in low gravity 
atmospheres. In the radiative zone it is unclear whether 
the vertical mixing is predominantly driven by convec- 
tive overshoot or gravity waves, but the picture emerg- 
ing from multi-dimensional hydrodynamical simulations 
in M dwarfs and brown dwarfs indicates that K zz is both 
dept h dependent and easily achieves values > 10 8 cm 2 
s" 1 (|Ludwig et al.ll2006t iFrevtag et al.ll2010t ). 

Figure [3] compares the best-fitting model photometry 
and spectrum to the observations. The model photom- 
etry agrees very well with the observations in nearly 
ever bandpass, with most bands agreeing at ler. The 
model J-band spectrum has about the right slope and 
nicely reproduces the water absorption starting at 1.33 
/im. At H band, the model also has roughly the correct 
shape but is slightly too linear across the central wave- 
lengths. The CO band in the K band is very well repro- 
duced by the model but the central wavelength region is 
aga in too flat. Skcmcr et al l (1201 It) also compare mod- 
els (|Burrows et alJl2006t iMadhusudhan et al.ll2011| ) with 
similar T c g and log (<?) to the same observations; how- 
ever, these models likely underestimate the non-LCE, as 
they do not adequately reproduce the near-IR spectrum, 
especially the CO band at 2.3 /im. An LCE version of 
the best-fit non-LCE model is shown in Figure [3] and 
demonstrates the impact non-LCE has on the near-IR 
spectrum. 

Also shown in Figure [3] is the 1600 K "Dusty" 
(jAllard et al.l 1200 1| ) equilibrium cloud model often se- 
lected as the best match to the near -IR spectra 
(jMohantv et al.l 120071: IFatience et al.ll2010ft . The 1600 
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Fig. 2. — Atmospheric properties for the best fitting model 
for 2M1207b. Top: temperature-pressure structure compared to 
condensation curves for two abundant cloud species. Middle: CO 
and CH4 mole fractions for equilibrium (dashed, red) and non- 
equilibrium (solid, with K z7j = 10 8 cm 2 s" 1 ) chemistry. Chemical 
and mixing timescales are also plotted (dotted lines). Bottom: 
dust-to-gas ratio for the intermediate cloud model (ICM) and the 
pure equilibrium cloud model. 



K model photometry does not agree well with the ob- 
servations across the full wavelength range. The 1600 K 
model spectrum, however, is very similar to the best-fit 
1000 K model across the K band, but only when the 
near-IR bands are scaled to match individually (to ac- 
count for the photometric disagreement for the 1600 K 
model). In both J and H bands, the 1600 K model over- 
predicts the peak flux and is noticeably more triangular 
at H than the observations. At K band, the 1600 K 
model provides only slightly better agreement with the 
observations than the 1000 K model. 

The remaining differences between the 1000 K model 
and near-IR spectra can be attributed to an incorrect 
proportion of dust opacity relative to molecular opac- 
ity. Given the simplicity of the cloud model used here 
such discrepancies are not surprising. Without a doubt, 
a more parameter-rich cloud model could be used to fine 
tune the comparison, but it is unlikely that such an ex- 
ercise will lead to significantly greater insights into the 
physical properties of the cloud. The primary lesson from 
this comparison is that atmospheric clouds and chemistry 
can dramatically alter the spectral shape and potentially 
lead to errors in effective temperature as great as 50%. 

The model comparisons to the photometry provide 
a new estimate for the bolometric luminosity. The 
mean luminosity, found by comparing to pure equilib- 
rium models, ICM/non-LCE models and black bodies, is 
logL/L Q = —4.68 with rms of 0.05. This luminosity is 
consistent with the earlier value mentioned above (based 
on if-band bolometric corrections). 

5. DISCUSSION AND CONCLUSIONS 

The best fitting T cS , log(#), and radius (1000 K, 10 4 cm 
2 , and 1.5 Ri U n) for 2M1207b are consistent with the 
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Fig. 3. — Photometric (top) and spectroscopic (bottom) comparison between 2M1207b observations (shown in black, see the text for 
details) and best-fitting model (blue). For comparison, synthetic photometry and spectra from a 1600 K model (red) with an equilibrium 
cloud model (aka "Dusty") along with an LCE model (green) with the same parameters as the non-LCE model. Surface gravities, radii, 
and effective temperatures are indicated in the legend. All fluxes have been scaled to 10 pc. 



demonstrates that including typical cloud thickness and 
non-LCE are all that is required to reproduce the current 
observations of 2M1207b. Such a model reminds us that 
the spectra of brown dwarfs are not strictly a function of 
temperature and, at young ages, can deviate significantly 
from expectations derived from older field brown dwarfs. 
The primary evidence supporting the edge-on-disk and 
protoplanetary collision hypotheses, was the previously 
deduced 1600 K effective temperature. A model with 
this temperature is shown to inadequately reproduce the 
available photometry and compares no better to near-IR 
spectra than a cooler cloudy model. Also, the disk-model 
comparison by ISkemer et al.l (|2011l ) further weakens the 
case for an edge-on disk. Consequently, no strong evi- 
dence remains for the previous disk or collision hypothe- 
ses. This conclusion is independent of the existence of 
HR8799b, but is certainly supported by it. 

The primary atmospheric contributors to the L-type 
appearance of 2M1207b, despite its low T e g, are clouds 
extending across the photosphere, thereby reddening the 
near-IR colors, and non-equilibrium chemistry, establish- 
ing a CO/CH4 ratio that is nearly the reciprocal of what 
is present in the photospheres of older field T dwarfs. The 
cloud properties, in particular the thickness, are proba- 
bly not substantially different from those found in late 
L dwarfs and the required K zz (~ 10 8 cm 2 s -1 ) is wel l 
within the range of field dwarfs. ISkemer et all (jMH) 
stress clouds as the primary explanation for the photo- 
metric and spectroscopic properties of 2M1207b. How- 
ever, non-LCE plays an equally important role. Proba- 
bly the most important underlying distinction between 
2M1207b and field brown dwarfs is low surface gravity, 
provided by its youth and low mass. Without low surface 
gravity, it is unlikely that clouds or non-LCE would be 
sufficient to give 2M1207b its current appearance. Such 
objects, therefore, should not be considered members of 
a new class, but rather represent the natural extension 



of substellar atmospheres to low gravity. 

Given their similar masses (~ 5Afj up ), it is possi- 
ble that 2M1207b and HR8799b represent two distinct 
states in the evolution of substellar atmospheric prop- 
erties (despite potentially different formation scenarios). 
After ~ 20 Myr of cooling, perhaps 2M1207b will spec- 
troscopically evolve into something resembling HR8799b 
and, eventually, into a traditional looking methane-rich 
T-dwarf. The very distinct spectra of the two objects 
(see Figure 15, Barman et al. 2011) would suggest rapid 
spectral evolution in the first 50 Myr, post formation. It 
is worth noting that the best-fit models for 2M1207b and 
HR8799b, though similar in atmospheric parameters, dif- 
fer in one significant respect - the radius derived from the 
effective temperature and luminosity. Our model here 
gives a radius of 1.5 -Rj UD for 2M1 207b, while even the 
coldest fit to HR8799b in iBarman etaD (|2011l ) gives a 
radius of only 1 i?j up . There are several possible expla- 
nations for this. First, the recovered radius is of course 
very sensitive to the effective temperature - a 100 K in- 
crease in 2M1207b and a 100 K decrease in HR8799b 
would remove the discrepancy, though neither would be 
a good fit to the spectra. Second, this may represent 
a fundamental difference in their internal state due to 
different formation scenarios. 2M1207b almost certainly 
did not form through a core-accretion process as it is 
highly unlikely that a disk surrounding a low-mass brown 
dwarf would have had enough material to accrete a giant 
planet, especially at large separations and in a short time. 
It is more likely that 2M1207b represents the tail end of 
binary star formation and, thus, might be expected to 
follow a similar cooling evolution as brown dwarfs. The 
formation of the HR8799 planets is less clear but could 
potentially have involved a n accretion per i od. T he "cold 
start" accretion models of Mar lev et al.l (|2007t ) predict 
significantly smaller radii for a given age and mass. Al- 
though the temperature and luminosity of HR8799b are 
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too high for the extreme cold start models, the smaller 
radius may be pointing toward a formation process that 
involved at least some loss of entropy. 

Finally, one can speculate on the implications that 
2M1207b and HR8799b might have on the spectral prop- 
erties for the broader young brown dwarf and planet 
population. If one adopts, ~ 1500 K as the upper T e g 
limit for T-dwarfs, then all T dwarfs younger than ~ 100 
Myr should be in the planet mass regime (< 13 Mj up ) 
and should have very low gravity (\og(g) < 4.5). How- 
ever, if 2M1207b and HR8799b are representative of cool, 
low gravity, substellar atmospheres, then non-LCE (and 
possibly clouds) will diminish the strength of CH4 ab- 
sorption across the H and K bands, making very young 
methane dwarfs rare. This prediction, however, is at 
odds w ith the tentative discoveries of ~ 1 Myr-old T 
dwarfs (IZapatero Osorio et alj|2002t iBurgess et alj|200& 
iMarsh et al.ll2010D . While at least one of these discover- 
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